Spin-split states in aromatic molecules and superconductors 
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A state where spin currents exist in the absence of external fields has recently been proposed to 
describe the superconducting state of metals. It is proposed here that such a state also describes 
the ground state of aromatic molecules. It is argued that this point of view provides a more natural 
explanation for the large diamagnetic susceptibilities and NMR shifts observed in these molecules 
than the conventional viewpoint, and it provides a unified description of aromatic molecules and 
superconductors as sought by F. London. A six-atom ring model is solved by exact diagonalization 
and parameters in the model where a ground state spin current exists are found. We suggest that 
this physics plays a key role in biological matter. 

PACS numbers: 



I. INTRODUCTION 



The large anisotropic diamagnetic susceptibility exhib- 
ited by aromatic molecules and the associated NMR fre- 
quency shifts at nearby atoms are generally interpreted 
as arising from delocalized ring currents in the molecu- 
lar rings [l|. This interpretation relies on molecular or- 
bital theory (Huckel theory) of tt electrons and assumes 
electrons are non- interacting 0, y]- On the other hand, 
calculations using valence bond theory taking electron- 
electron interaction into account predict that the tt elec- 
trons in aromatic rings are localized [J]. These electrons 
move in a one-dimensional half-filled band, where the 
effect of electron-electron interactions should be strong; 
electron-electron interactions in a one-dimensional half- 
filled band can lead to localization of electrons (Mott 
insulating state) which would suppress the diamagnetic 
response. Within valence bond theory the aromatic char- 
acter is proposed to result from the coupling of electron 
spins |J|. 

F. London suggested long ago that the diamagnetic 
currents in aromatic rings are analogous to supercurrents 
in superconductors 5j. At that time superconductivity 
in metals was not understood from a microscopic point 
of view. After the advent of BCS theory the possible 
connection between ring currents in aromatic molecules 
and superconducting currents has been discussed by W. 
Little [6] in connection with a proposed excitonic pairing 
mechanism. However, to date no compound has been 
identified that would exhibit Little's excitonic supercon- 
ductivity. 

The present author suggested in 1990 that the large 
diamagnetic response of aromatic ring molecules may 
originate in the existence of a spin current in the ground 
state of such molecules in the absence of applied fields [3] ■ 
Thus, an applied magnetic field would not set electrons 
into motion but rather slow down the component of the 
pre-existent spin current in one direction and speed up 
the spin current in the opposite direction. A heuristic 
model for this physics was proposed and it was suggested 
that it provides a natural explanation for various exper- 
imental observations. At that time we did not make a 



connection between this physics and superconductivity. 
However our recent work on the theory of hole 
superconductivity's] led to the conclusion that metals ex- 
pel negative charge from the interior towards the surface 
in the transition to the superconducting state [9]. This 
phenomenon (which also provides a 'dynamical' expla- 
nation of the Meissner effect [io| ) leads to the prediction 
that a a macroscopic spin current exists near the sur- 
face of superconductors in the absence of applied exter- 
nal fields (Spin Meissner effect) [ill. [l3|. Naturally this 
raises the question whether the spin currents in aromatic 
molecules proposed in ref.Q are related to the spin cur- 
rents in superconductors recently predicted [l2|. In this 
paper we propose that in fact such a deep connection 
exists, validating the early intuition of F. London. A 
state qualitatively different from the Huckel[2,3 and va- 
lence bondQl states, carrying a spin current around the 
ring, is proposed to describe the ground state of aromatic 
molecules. 



II. SPIN CURRENT AND tt FLUX IN 
SUPERCONDUCTORS 

We have recently proposed that the Meissner effect in 
superconductors can be understood 'dynamically' if elec- 
trons expand their orbits to mesoscopic orbits of radius 
2Al in the transition to superconductivity, with Ai the 
London penetration depth [l^. In the absence of applied 
magnetic field this orbit expansion results, through the 
spin-orbit interaction of the electron magnetic moment 
/is — eh/2meC with the positive ionic background \ e\ns 
{ug =superfluid density), in a spin current with speed [12| 
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and opposite direction for electrons with opposite spin. 

The angular momentum of electrons in orbits of radius 
2Al with speed Eq. (1) is 



L = mef°(2Ai) = 



(2) 



This remarkable result indicates that electrons in super- 
conductors have an orbital 'spin' analogous to the intrin- 
sic electron spin. A spinning electron can be visualized 
as a charge orbiting at speed c in a circle of quantum 
electron radius 
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and the electron in the superconductor can be interpreted 
as an 'amplified image' of the spinning electron, with 
amplification factor 2Af,/r^[f3l|. 

An angular momentum h/2 corresponds to a wavefunc- 
tion ^ = |^|e*^ depending on the azimuthal angle (p as 



* oc e'^^f^ = e^*'^/2 



(4) 



so that the electron acquires a phase tt in traversing a full 
orbit. This in turn can be understood as arising from a 
'flux' 00 inside the orbit giving rise to the phase factor 
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with 00 = tf ^^^ ^^^^ quantum and Ag. the spin-orbit 
vector potential [ij]. 

Thus, the development of the macroscopic spin cur- 
rent predicted by the theory of hole superconductivity 
upon the establishment of macroscopic phase coherence 
in superconductors requires 'closing' [15[ of the wavefunc- 
tion for each member of the Cooper pair according to the 
relation 



(j) V0a ■dl = ±Tr 



(6) 



where 9^ is the phase of the member of the Cooper pair of 
spin a, instead of the usual Bohr-Sommerfeld quantiza- 
tion rule. Alternatively we can think of 9^ as character- 
izing the phase of the entire condensate of spin a rather 
than a single member of a Cooper pair, since the phases 
of different Cooper pairs are phase locked due to macro- 
scopic phase coherence. The phase factor e**'^ = — 1 
acquired by one spin component is cancelled by the same 
factor in the opposite spin component so that the super- 
conducting wavefunction is single- valued. 

As is well known, for an individual spin-1/2 spinor, 
a rotation by 360 degrees gives rise to a change in 
sign 16], which we can also interpret as arising from an 
enclosed tt— flux. Thus, a single spinning electron shows 
the same phase behavior as the Cooper pair member in 
the mesoscopic orbits in the superconductor within the 
theory of hole superconductivity, proposed to describe 
all superconductors [17]. The fact that the same physics 
shows up at the subatomic scale Vq (Eq. (3)) and at the 
mesoscopic scale 2Xl suggests that it will also show up at 
other length scales provided phase coherence exists. In 
this paper we propose that this physics shows up at the 
atomic scale ao = h^ /rriee^ ^ ^rg(2Ai) and in particu- 
lar manifests itself in the existence of a spin current in 
the ground state of aromatic molecules such as benzene. 







(a) 



(b) 



FIG. 1: The two possible orientations of the carbon p^ orbitals 
in the benzene molecule. In the conventional point of view 
only one of them (either one) plays a role. 




FIG. 2: In the point of view proposed here, as the electron 
goes around the benzene ring once, it ends up in the orbital 
with opposite sign. Two rounds are needed to go back to the 
original state. 



III. SPIN CURRENT AND n FLUX IN 
BENZENE 

The aromatic character of benzene is ascribed to the 
six electrons in the pn orbitals oriented perpendicular to 
the plane of the molecule. There are two possible orien- 
tations for the atomic wave function, as shown schemat- 
ically in Fig. 1. In the conventional point of view, these 
two orientations are equivalent except for a sign conven- 
tion which has no physical significance, and either one 
of them is used. Instead, we propose here that the fact 
that there are two possible orientations of the Pz orbital 
as shown in Fig. 1 has profound physical significance. 

We propose, by analogy with the situation for the in- 
trinsic electron spin 16], as well as for the 2Al orbits 
in the superconductor discussed in the previous section, 
that as an electron goes around once in the aromatic ring, 
it ends up in the orbital with opposite sign orientation, 
as shown schematically in Figure 2. That is, the elec- 
tron acquires a tt phase shift, and two rounds are needed 
to bring the electron back to its original state with the 
same sign. Because for every spin-up electron there is a 
spin-down electron moving around in opposite direction, 
no sign change occurs for the wavefunction of the system 
as a whole, just like in the case of the superconductor. 

To describe this scenario we may take the tight binding 
Hamiltonian kinetic energy with hopping between near- 
est neighbors < ij > 
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and take hoppings Uj — t for all bonds except one, where 
tij — — t, corresponding to 'antiperiodic' boundary con- 
ditions. Thus, the electrons acquire a phase tt in going 
around the ring. The non-interacting spectrum and the 




the spinor[l| 



FIG. 3: In the conventional scenario (a), the ground state of 
benzene is non-degenerate and there is no spin current. In 
the scenario proposed here (b) the (non-interacting) ground 
state is four- fold degenerate and can carry spin current. 




FIG. 4: In the scenario proposed here, the non-interacting 
ground state is four-fold degenerate as shown in the Figure. 
Two of the ground states, (a) and (b) , carry spin current and 
the other two, (c) and (d), carry charge current. 



occupation of the states for 6 electrons is shown in Fig. 
3 contrasted with the conventional case. In the scenario 
considered here there are four degenerate ground states 
in the non-interacting case, shown in Fig. 4. Two of 
them carry a spin current. 

The energy level structure shown in Fig. 4 would arise 
within the conventional viewpoint if there is an applied 
magnetic field with flux (Jjq within the aromatic ring. 
Here we propose that this structure arises spontaneously 
in benzene in the absence of applied magnetic field. It 
can be thought of as arising from n — /lux of opposite 
sign affecting the states of up and down electrons and 
parallels the behavior predicted for superconductors by 
the theory of hole superconductivity [12]. 

Alternatively we argue as follows. Consider a rotation 
of the coordinate system by an angle ip — n/S around the 
0— axis chosen perpendicular to the plane of the aromatic 
ring. This is a symmetry operation that leaves the ben- 
zene ring invariant. A spinor x localized at a given atom 
in the ring is described in the new coordinate system by 



X 
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with 9 — 7r/6. Six such transformations in sequence bring 
the coordinate system back to the original one, however 
the spinor changes its sign. We now consider these sym- 
metry operations from an "active" rather than a "pas- 
sive" point of view. Namely, instead of rotating the ob- 
server's coordinate system we rotate our physical system, 
the spinor. If x is an eigenstate of a^ , it acquires a phase 
factor e^''^/^ in hopping from a site to a neighboring site, 
with opposite sign phase factor for both eigenstates of (7^ ■ 
Adopting this point of view, the kinetic energy opera- 
tor is, instead of Eq. (7) 
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(9a) 



(9b) 



(9c) 



and the spin quantized in direction perpendicular to the 
plane of the molecule. The energy spectrum is identical 
to that shown in Figure 4. In Eq. (9a), < ij > indicates 
nearest neighbor atoms with i to the left of j. 



IV. ELECTRON-ELECTRON INTERACTIONS 

Can states that carry a spin current, such as those 
shown in Fig. 4 (a) and (b), be stable in the presence 
of strong electron-electron interactions? It may be sus- 
pected that in a finite ring electron-electron scattering 
will destroy the spin current by coupling states with spin 
current flowing in opposite directions |18j. 

We consider the kinetic energies Eq. (7) or (9) in the 
presence of on-site and nearest neighbor Coulomb repul- 
sion as well as nearest neighbor ferromagnetic exchange 
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+ E jf4aCj.cl^-C..' (10) 

<i.j>aa' 

With the kinetic energy Eq. (7) we have simply Jf,"^ = 
J > 0. This term (ferromagnetic nearest neighbor ex- 
change term) arises from an off-diagonal matrix element 
of the Coulomb interaction in the tight binding formula- 
tion (J) which is always positivc[19]. It was proposed 
in Ref.[7| that it plays an essential role in stabilizing 
the spin-split state in aromatic molecules. However, in 
Ref. lZI , the spin current state envisioned was based on ex- 
citations of the conventional non-interacting system Fig. 
3(a). We find here that J also plays an essential role in 
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FIG. 5: Spin current versus nearest neighbor ferromagnetic 
exchange interaction strength for three values of the staggered 
flux 4>s in Eq. (12). 



stabilizing the spin current state based on the scenario 
of Fig. 4. If we use the kinetic energy Eq. (9), the J 
interaction in Eq. (10) takes the form 






(11) 



so that it acquires a phase factor for exchange of electrons 
of opposite spin. The spectrum of the Hamiltonian Eq. 
(10) is the same for the two forms of the kinetic energy 
Eq. (7) and Eq. (9) and we will use Eq. (9) in what 
follows. 

We diagonalize the Hamiltonian Eq. (10) exactly for 6 
electrons with equal number of up and down spins. It is 
easy to see that in the presence of interactions the ground 
state cannot carry a spin current if it is non-degenerate, 
since such a spontaneously broken symmetry state cannot 
exist in a finite system. So we look for parameters in the 
Hamiltonian that will give rise to a doubly degenerate 
or nearly doubly degenerate ground state. Indeed we 
find that generically such Hamiltonian parameters exist 
even for the case of strong interactions. We find them in 
the neighborhood oi U ^ 2V and they require non-zero 
values of the exchange interaction J. 

Consider for example the case U = 4, V = 2 (with 
t = 1). To detect a spin current we add to the phases of 
the spinors a small staggered fiux (f>s, i.e. take 



= {-+(j)s)<J. 

6 



(12) 



This arises in the real system from the spin-orbit inter- 
action not included in the Hamiltonian Eq. (10), as dis- 
cussed in the next section. Figure 5 shows the ground 
state spin current 
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FIG. 6: Four lowest energy levels of the Hamiltonian for 
U = A, V = 2 as function of J. The ground state switches 
from singlet (s) to triplet (i) at J ~ 0.714. At that point 
(indicated by the symbol in the figure) the degenerate ground 
state corresponds to the spin-current carrying states of Fig. 
4 (a), (b). 
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FIG. 7: Four lowest energy levels of the Hamiltonian for 
(7 = 8, y = 3.5 as function of J. The singlet-triplet crossing 
point where the spin current is largest is at J ~ 0.361. 



versus the nearest neighbor exchange amplitude J for 
three small values of the staggered flux (/)„. It can be 
seen that a spin current exists for finite staggered fiux 
and non-zero J . As the staggered flux approaches zero, 
the spin current approaches a (5— function at one value 
of J, J = 0.714 for this case. This corresponds to the 
crossing of the two lowest energy levels, a singlet and a 
triplet, as shown in Fig. 6. In the limit (ps -^ 0, spin 
current will only exists at the precise point where the 
energy levels cross. 

Similarly we show in Figs. 7 and 8 the four low- 
est energy levels and the ground state spin current for 
symmetry-breaking field (pa = 0.0001 for two other sets 
of [/, V values. It can be seen that the behavior is generic. 
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FIG. 8: Four lowest energy levels of the Hamiltonian for 
U — 2, V = 1 as function of J. The singlet-triplet crossing 
point where the spin current is largest is at J ~ 0.192. 
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FIG. 9: Four lowest energy levels of the Hamiltonian for 
U = 4, V = 1.85 as function of J. The singlet-triplet crossing 
point where the spin current is largest is at J ~ 0.293. 



Note that the range of parameters in the model where 
we find this interesting behavior is quite reasonable. For 
example, first principles calculations for benzene yield 
values t = 2A0eV, U = 11.26ey;2l]. The nearest neigh- 
bor repulsion V estimated from the Ohno formula[22l| 
V{R) = U/y/TTWiPje^ yields V == 7.6eV^ for near- 
est neighbor distance R = lAA. Hence U/t = 4.7, 
V/t = 3.2. The estimation of these parameters is subject 
to ambiguity, for example related to the validity of the 
zero differential overlap approximation [2^. Furthermore 
the parameter values could be renormalized by inclusion 
of further neighbor interactions in the model. We will 
investigate these questions in separate work. 

The range of nearest neighbor interaction values where 
this behavior occurs is rather restricted. Fig. 9 shows the 
lowest energy levels ioi U = A, V = 1.85. The crossing 
point where spin current exists for infinitesimal spin-orbit 



FIG. 10: Four lowest energy levels of the Hamiltonian for 
U = 4, V = 1.82 as function of J. The singlet-triplet crossing 
point is at J ~ 0.143 indicated by the symbol in the figure, 
however another singlet energy level lower in energy exists for 
that J that does not carry a spin current. 
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FIG. 11: Charge current response to applied magnetic field 
relative to non-interacting case with periodic boundary con- 
dition, for the case (7 = 4, V = 2. The conventional scenario 
(no flux) gives somewhat larger charge current than the sce- 
nario proposed here. The small kink in the tt— flux case near 
J = 0.7 corresponds to the energy level crossing point in Fig. 
6. 



coupling has moved from J = 0.714 down to J = 0.293. 
For smaller V the crossing point for these energy levels 
moves further to the left and crosses another energy level, 
that becomes the lowest energy state and does not carry 
spin current, as shown in Fig. 10 for U = A, V = 1.82. 
However even in this situation if the spin-orbit interac- 
tion is finite the ground state and first excited state will 
carry a spin current for J larger than the crossing point 
with the non-spin-current carrying ground state (J > 0.2 
in Fig. 10). 

Thus, for U — A the range of V where a spin cur- 
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FIG. 12: Electrons move in an outward-pointing electric field 
E originating in the charge transfer from the C atom to the 
C — H bond. The spin down (up) electron has lower energy 
when orbiting counterclockwise (clockwise), as shown in the 
figure, due to spin-orbit coupling. 



rent carrying ground state for an infinitesiinal </><, exists 
is restricted to 1.83 < V < 2.05. Similarly restricted 
parameter ranges where spin current in the ground state 
can develop for infinitesinial spin-orbit coupling exist for 
other U values. 

The response of the system to an applied external mag- 
netic field is similar to the conventional scenario for the 
case considered here. Figure 11 shows the charge current 
that develops for the case [/ = 4, F = 2 as function of 
J, relative to the current expected within Huckel the- 
ory (non-interacting electrons) in the presence of a small 
magnetic flux, of the same sign for both spin orienta- 
tions. Both the conventional scenario with no spin-orbit 
flux and the scenario with tt— flux proposed here give a 
charge current, hence a diamagnetic susceptibility, that 
is almost constant as function of the interaction J and 
not much smaller than the non-interacting value (which 
is 1 in Fig. 11). The small glitch in the lower curve 
corresponds to the point where the energy level crossing 
occurs in Fig. 6. 



V. SPIN-ORBIT COUPLING 

If the ground state of the interacting Hamiltonian is 
degenerate (or nearly degenerate) in the way discussed 
in the previous section, spin-orbit coupling will split the 
degeneracy and give rise to two spin-current-carrying 
low-lying states. The lowest of the two states is shown 
schematically in Fig. 12. 

We can understand the physics as follows. The electron 
in the spin current moves in an electric field pointing out- 
wards, just like in an atom or in the superconductor ac- 
cording to the theory of hole superconductivity [12t|. The 
spin-orbit interaction Hamiltonian in electric field E is 
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and it gives rise to an energy splitting 

eh 



AE 



2mic- 



■Ep 



(15) 



(16) 




FIG. 13: Allowed optical transitions in benzene in the con- 
ventional scenario ((a)) and in the scenario proposed here 
((b)) are indicated by the dashed lines. 



between electrons with spin parallel and antiparallel to 
Exp. Taking p ^ h/a, with a the radius of the orbit (or 
equivalently the C — C distance) 
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eh^ 



2m'lc^a 
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(17) 



The spin-orbit interaction gives rise to the Aharonov- 
Casher vector potential [23l| 



2e 



(18) 



with fl — /is (7, fiB the Bohr magneton. In going around 
the ring the electron acquires a phase shift 



seAc = ^ 

he 



A„ 



„ 27rae 
dl = -^ — A„ = 



he " 2TOeC^ 

or in terms of the energy splitting Eq. (17) 
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Assuming the energy splitting AE is given by the atomic 
spin-orbit splitting in C between ^P3/2 and ^Pi/2 states, 
AE = 8meV, gives for the AC phase Eq. (20) SOac = 
0.065. The staggered flux introduced in Eq. (12) is then 
l/6th of this value, i.e. 



= 0.0011 



(21) 



close to the middle value oi (f)s used for illustration in Fig. 
5. For spin-orbit energy splitting a factor of 10 larger, 
i.e. AE ^ SOmeV, the flux corresponds to the largest 
value shown in Fig. 5, where a rather large spin current 
occurs in a wide range of J values. 

In summary, we propose that the two lowest energy 
states of benzene carry a large spin current around the 
ring, as depicted schematically in Fig. 4 (a) and (b), 
and have a small energy separation of order meV . The 
ground state has spin current direction as shown in Fig. 
12. 



VI. PHYSICAL CONSEQUENCES 

One immediate consequence of the scenario proposed 
here concerns optical properties and is illustrated in Fig. 



13. As discussed by Piatt [24|, allowed optical transitions 
in benzene require change of one unit of angular momen- 
tum, corresponding to change of the electron wave vector 
by 27r/6. In the conventional scenario there is only one 
way to do this starting from the ground state, as shown 
schematically in Fig. 13 (a). In the scenario proposed 
here there are two distinct ways of exciting electrons to 
states differing in wave vector by 27r/6, as shown in Fig. 
13 (b). The photon wavelengths associated with these 
two processes should be different because of the different 
effects of electron-electron interactions in the final states. 

Experimcntally[25] three bands are seen in the op- 
tical absorption spectrum of benzene, at wavelengths 
around A = 2600A, 2050A and 1850A. The one at 
2600 is weak and has considerable fine structure, thus 
it is likely to arise from 'forbidden transitions' involv- 
ing interaction with vibrational degrees of freedom |26l|. 
The other two bands are much stronger, and their inter- 
pretation has been controversial, because there is only 
one allowed transition in the conventional scenario (Fig. 
13(a)). There is general agreement that the transition at 
A = 1850^ is an allowed one because it is the strongest, 
and the transition at A = 2050A has been attributed to 
one or another forbidden transition. A variety of calcu- 
lations of oscillator strength have been reported [27H30J 
but it is not clear that these complicated calculations 
explain the experimental observations. Instead, our sce- 
nario predicts two allowed transitions (Fig. 13(b)) and 
hence two strong absorption bands in the optical absorp- 
tion of benzene, thus offering the possibility of a more 
natural explanation for the observations. 

The ground state and low-lying energy state proposed 
here for benzene are neither pure singlet nor triplet, in- 
stead singlet and triplet are strongly mixed, with nearly 
equal weights if J is near the crossing point of singlet and 
triplet energy levels shown in the figures. It is possible 
that transitions from higher excited states to these states 
will provide a more natural explanation of the phospho- 
rescence properies of benzene than the conventional un- 
derstanding, within which phosphorescence in benzene is 
'doubly forbidden' and requires vibronic interactions to 
play a key role|3l|. 

The spin current proposed to exist in benzene will give 
rise to an electric field [3]. A moving magnetic moment 
jl is equivalent to an electric dipole moment p in the 
laboratory frame 



p= - X fj. 
c 



(22) 



and will give rise to an 'effective' dipolar charge distri- 
bution depicted schematically in Fig. 14. This dipolar 
charge distribution generates a quadrupolar electric field 
around the aromatic ring. For a point in the plane at 
distance x » a from the center of the ring 



E{x) ^ — 



(23) 





a-state 



p-state 



with a the radius of the ring. At distance a: = 2a of the 
center, the electric field \s E ^ p/o?- From Eq. (22) 



FIG. 14: Dipolar charge distribution produced by spin cur- 
rent. Q— state: spin antiparallel to orbital angular momen- 
tum, /?— state: spin parallel to orbital angular momentum. 
I/CT denotes the orbital angular momentum for spin a and only 
the upward pointing La is shown. The q— state will generally 
be lower in energy due to spin-orbit coupling. 



the dipole moment is p '^ 77^VA^ for benzene, and the 
magnitude of the electric field aX x = 2a is E ^ bAfiV/A. 

We suggest that the two low-lying states of the aro- 
matic ring depicted in Fig. 14, which we call a and (3 
states, are stable, because of the topological constraint 
that the 'winding number' for electrons going around the 
ring is 1/2, and that a large energy barrier exists for 
switching between a and /3 states. If these states exist 
in aromatic rings they are likely to play a crucial role in 
biological matter where aromatic rings are ubiquitous. 

How can an aromatic ring 'switch' from one to the 
other low-lying state? The /3 state becomes the lower en- 
ergy state if the electric field that the electrons in the 
spin current 'see' points inward rather than outward, 
which could occur upon changes in the electrostatic field 
around the ring caused by chemical reactions involving 
charge transfer, or transport of ions or electrons. A pro- 
ton going through the interior of a benzene ring [32] may 
switch a /3 to an a configuration due to the additional 
strong electric field pointing outward, and conversely an 
electron going through the center of a benzene ring may 
switch the a to the /3 state. 

It is tempting to speculate that the a and /3 states 
could act as quantum qubits in a topological quantum 
computer [33J, namely the brain. Electric fields origi- 
nating in the spin current, of the magnitude estimated 
above, and resulting electric potential differences may be 
involved in neural signal transmission. Or, the a and /3 
states in aromatic rings in the brain may just be stor- 
ing one classical bit of information and play a key role 
in memory storage [7|. Phase coherence of spin ring cur- 
rents over large distances could play an important role in 
various aspects of brain activity such as sleep. More gen- 
erally, these stable dynamic low-lying states of aromatic 
rings may play a key role in a wide variety of biological 
processes in living organisms [34J]. 

As discussed earlier, the theory of hole superconductiv- 
ity predicts that superconductors in the ground state also 
carry a spin current and associated tt flux. The ground 
state of a small superconducting ring corresponds to the 
a state in Fig. 14. It should be possible to switch such a 
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ring to the (3 state, perhaps through transient apphcation 
of a large electrostatic field or a very large magnetic field. 
Such ring would exhibit a stable persistent spin current 
in opposite direction to that of the ground state (in the 
absence of external fields), and experimental detection of 
such a state would provide convincing proof of the the- 
ory discussed here and in ref . fl2] . The existence of such 
spin currents in superconducting rings could be experi- 
mentally demonstrated by detection of the quadrupolar 
electric fields generated by the dipolar charge distribu- 
tions shown in Fig. 14, or by measuring the resulting 
force between small superconducting rings with spin cur- 
rents, as discussed in ref. 12|. 



VII. DISCUSSION 

Aromatic rings have An + 2 atoms, with n integer. 
It is generally believed that such molecules have a sin- 
glet ground state described by delocalized electrons in 
a half-filled band [3, Q, with a large energy gap to the 
first excited state. Alternatively, it has been proposed 
that a similar energy level structure and aromatic char- 
acter can be described by a valence bond state with lo- 
calized electrons and coupling between electron spins [J|. 
In this paper we have proposed instead a qualitatively 
different scenario: that the ground state of aromatic ring 
molecules is a mixture of singlet and triplet states that 
carries a spin current around the ring, and that another 
such state exists nearby in energy carrying spin current 
in the opposite direction. We have also proposed that 
the same is true for small superconducting rings made of 
ordinary superconductors. 

This scenario follows from the fundamental assump- 
tion that when an electron goes around a closed path in 
a coherent way it acquires a tt phase shift, and two rounds 
are needed to bring it back to the original state, like a 



point in a Moebius strip. In the aromatic ring, starting 
in a Pz orbital the electron ends up in the (—Pz) orbital 
after one round. We have justified this assumption by 
invoking the transformation properties of the spinor un- 
der rotation in quantum mechanics [16|, and by analogy 
with the phase behavior predicted by the theory of hole 
superconductivity, stemming from the finding that elec- 
trons in a Cooper pair carry orbital angular momentum 

Within this assumption, and using a Hubbard-like 
Hamiltonian including nearest neighbor exchange to 
model the benzene ring, we have found reasonable pa- 
rameter ranges in the interacting system where a spin- 
current-carrying ground state and low-lying excited state 
exists. The same is likely to be true for other aromatic 
ring structures using the same model. We have argued 
that our scenario provides a more natural explanation for 
the optical properties of benzene as well as for the exis- 
tence of ring currents evidenced by its large diamagnetic 
susceptibility and proton NMR shifts. 

Whether our proposal is true or false can be decided 
experimentally. It should be possible to make aromatic 
rings and small superconducting rings in both predicted 
low-lying spin current states and study their properties. 

If our proposal is true it is likely to yield fundamen- 
tal insights into the properties of biological matter, as 
well as on quantum mechanics itself, changing the cur- 
rently accepted understanding in fundamental ways. It 
is also likely to open up possibilities of qualitatively new 
technological advances such as quantum computation. 



Acknowledgments 

The author is grateful to Congjun Wu for stimulating 
discussions. 



J. Wiley, New York, 1986 



[1] P.J. Garratt, 'Aromaticity' 

and references therein. 
[2] L. Pauling, J. Chem. Phys. 4, 673 (1936). 
[3] R.C. Haddon, J.Am.Chem. Soc. 101, 1722 (1979). 
[4] D.L. Cooper, J. Garratt and M. Raimondi, Nature 323, 

699 (1986). 
[5] F. London, J. Phys. Radium 8, 397 (1937). 
[6] W. A. Little, Phys. Rev. A134, 1416 (1964). 
[7] J.E. Hirsch, Mod. Phys. Lett. B4, 739 (1990). 
[8] J.E. Hirsch, J. Phys. Chem. Solids 67, 21 (2006) and 

references therein. 
[9] J.E. Hirsch, Phys.Rev.B 68, 184502 (2003). 
[10] J.E. Hirsch, J. Phys. Cond. Matt. 20, 235233 (2008); 
Jour, of Superconductivity and Novel Magnetism 23, 309 
(2010) 
[11] J.E. Hirsch, Phys.Rev.B 71, 184521 (2005) 
[12] J.E. Hirsch, Europhys. Lett. 81, 67003 (2008). 
[13] J.E. Hirsch, arXiv 1002.2688 (2010). 
[14] J.E. Hirsch, Ann. Phys. (Berlin) 17, 380 (2008). 



[is; 

[16 

[17 
[18 

[19 
[20 

[21 

[22' 
[23 

[24 
[25 

[26 

[27 



A.V. Nikulov , Phys. Rev. B 64, 012505 (2001); arXiv 

0812.4118 (2008); arXiv 1004.3914 (2010). 

E. Merzbacher, 'Quantum Mechanics', Third Edition, 

John Wiley and Sons, Inc, New York 1998, Chpt. 16. 

J.E. Hirsch, Physica Scripta 80, 035702 (2009). 

C.J. Wu, private communication. 

J.E. Hirsch, Phys.Rev.B 40, 2354 (1989). 

J.E. Hirsch, J.Chem.Phys. 171, 161 (1990). 

Z. G. Soos, Y.A. Pati and S. K. Pati, J.Chem.Phys, 112, 

3133 (2000). 

K. Ohno Theoret. Chim. Acta 2, 219 (1964). 

Y. Aharonov and A. Casher, Phys. Rev. Lett. 53, 319 

(1984). 

J.R. Piatt, J. Chem. Phys. 17, 484 (1949). 

J.N. Murrell, "The Theory of the Electronic Spectra of 

Organic Molecules", Wiley, New York, 1963. 

C. Herzberg and E. Teller, Z. Physik. Chem. B21, 410 

(1933). 

J.N.MurreU and J.A. Pople, Proc. Phys. Soc. A69, 245 



(1956). 1147 (1973). 

[28] T.M. Dunn and C.K. Ingold, Nature 176, 65 (1955). [32] R.S. Shresth et al, Curr. Sci. 71, 49 (1996). 

[29] A.C. Albrecht, J. Chem. Phys. 33, 156 (1960). [33] C. Nayak et al, Rev.Mod.Phys. 80, 1083 (2008). 

[30] P.M. Johnson, H. Xu and T.J. Sears, J. Chem. Phys. [34] P. Ehrenfreund, S. Rasmussen, J. Cleaves and L. Chen, 

125, 164330 (2006). Astrobiology 6, 490 (2006). 

[31] J. van Egmond and J.H. van dor Waals, Molcc. Phys. 26, 



